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Introduction 

^ , Long-lasting spiral structure in galaxies is established when density waves 

^ I form a pattern that rotates at a constant angular frequency Q,p, the pattern speed. 

■ This quantity, critical to understanding spiral structure, galaxy evolution, and 

^ I the connection, if any, between bars and spirals in galaxies, is not directly 

measurable. Most determinations of pattern speeds have relied on identifying 
Q . morphological features at certain radii as resonances, but this depends signif- 

^ I icantly on knowledge of the behavior of stars and gas at these resonances. 

Q^' The Tremaine-Weinberg (TW) method allows for the determination of pattern 

Q ! speeds based on kinematic observations in a way completely independent of 

V-j ' resonance identification, or theories of spiral structure, relying on the identi- 

c/3 . fication of a tracer that obeys the continuity equation as it orbits through the 

^ ' pattern. 

This work expands on previous studies of the TW method applied to galax- 
^ , ies with molecular- or Hl-dominated ISMs. In galaxies where the ISM is found 

Jh I to be dominated by one phase, it is assumed that this component satisfies the 

continuity requirement; for galaxies with H2 dominated ISMs, for instance, 
the TW method has successfully been applied using ^^CO(l-O) emission as a 
tracer for the molecular gas surface density (Rand & Wallin 2004). 

1. The Tremaine-Weinberg Method 

If the component obeys continuity as it orbits in the spiral/bar potential, 
Tremaine & Weinberg (1984) show that 
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where {X, Y) are the sky coordinates || and _L to the major axis; Y) is 
the surface density of the component; V|| is the line-of-sight velocity; and i is 
the incHnation angle of the galaxy. 

If the intensity of the component's emission traces the surface density, then, 
at a given Y (cf. Merrifield & Kuijken 1995), 

oo oo 
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and modulo sin i, can be determined from total intensity and velocity field 
maps by measuring intensity-weighted averages of observed position <x> and 
velocity <v> along lines of constant Y (parallel to the major axis). The slope 
in a plot of <v> vs. <x> for such "apertures" is the pattern speed. 

2. Data and Analysis 

Presently, to determine the spiral density wave pattern speed, we use high- 
resolution CO data cubes from the BIMA Survey of Nearby Galaxies (SONG). 
The BIMA SONG observations include single-dish data, thereby sampling all 
spatial frequencies to support the continuity requirement. Pattern speed de- 
terminations for three of the galaxies in this work are presented based on CO 
emission data cubes alone; we cite the results of Wong & Blitz (2002) wherein 
radial profiles of NGC 5055, 5033, and 5457 (MlOl) indicate that these galax- 
ies are molecule-dominated within the region of CO detection. 

In this work, a constant conversion factor X=2xl0^^ mol cm~^(K km 
s~^)~^ between CO intensity and H2 column density is assumed. As demon- 
strated in Zimmer, Rand & McGraw (2004), the choice of conversion factor 
has little effect on the TW calculation. The choice, however, clearly influences 
assertions of molecular dominance, particularly with regard to an important 
suspected variation of X with metaUicity, which is supported by empirical ev- 
idence for an inverse correlation (e.g. Boselli, Lequeux, & Gavazzi 2003). For 
the galaxies studied here, NGC 5457, 5055 and 5033, metaUicity estimates are 
subsolar (as compiled by Pilyugin, et al. 2004), indicating that, if anything, the 
conversion factor underestimates the amount of H2 present. 

IncHnation and position angles (PA's) used in the TW calculation are deter- 
mined by fitting tilted rings to the velocity field as performed by the GIPSY 
task ROTCUR. As found in Debattista (2003) and Rand & Wallin (2004), the 
scatter in the slope of the <v> vs. <x> plot worsens for PA's further from the 
nominal value. Error bars reported for the pattern speed of each galaxy account 
for PA uncertainty and the scatter in the <v> vs. <x> plot for nominal PA. 
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Morphological 


D 


i 


PA 


^sys 






Type 


(Mpc) 
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n 


(km/s) 


(km/s/kpc) 


NGC 5457 


SAB(rs)cd 


7.4 


21 


42 +4/-2 


257 


48 +7/-6 


NGC 5055 


SA(rs)bc 


7.2 


57 


99 ±1 


506 


67 ±6 


NGC 5033 


SA(s)c 


11.8 


68 


352 ±2 


880 


93 +18/-39 



Table 1. Properties of sample galaxies. 



3. Results 

Given that NGC 5457, 5055 and 5033 appear to be dominated by molecular 
gas in their inner disks, we can be reasonably confident that the continuity re- 
quirement of the TW method is satisfied. Table 1 lists the properties and results 
of the TW calculation for each of the galaxies in our sample. Figures la and 
b display an example of the line-of-sight apertures used in the TW calculation 
for NGC 5457 and the corresponding <v> vs. <x> plot, respectively. 

4. Conclusions and Future Work 

The Tremaine-Weinberg method has long been apphed with much success 
to bars in galaxies. The relatively recent applications of the method to spiral 
patterns, too, are compelling: pattern speeds for M51 and M83, for example 
(Rand & Wallin 2004), were found to be in good agreement with those de- 
termined from resonances. Generally, however, it is unclear how reliable spi- 




Figure la. First moment map of the 
BIMA CO data for NGC 5457 with aper- 
tures used in the TW calculation spaced at 
the resolution of the map shown overlayed 
as solid lines, with the westernmost being 
aperture 1. 



Figure lb. Plot of <v> vs. <x> for the 21 
apertures shown in fig. la. The dashed line 
is the best-fit straight line to all apertures. 
The corresponding value of fip is 48 +7/-6 
km s~^ kpc~^. 
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aperture # 



Figure 2 Plot of <v>/<x> 
vs. aperture number for NGC 
5457 with filled and open cir- 
cles representing values from 
the TW calculation using the 
original and axisymmetric 
velocity fields, respectively. 



ral pattern speed determinations with the TW method are, considering that the 
converging and diverging flows responsible for the pattern are typically weaker 
in a spiral (e.g. Roberts & Stewart 1987) than in a bar; the method, based on 
the continuity equation, reUes on the detection of such non-axisymmetric, or 
"streaming", motions that trace spiral arm surface density enhancements in 
order to calculate a pattern speed. It is easy to show, in the case of purely 
axisymmetric azimuthal motion, that the method returns a weighted average 
of the angular rotation frequency. This makes sense: in such a case, only a 
material pattern can be sustained. 

To understand whether streaming motions are significantly contributing to 
the TW analysis, we have repeated the procedure using axisymmetric veloc- 
ity fields built from the output of ROTCUR. We generally find Uttle difference 
relative to the result for the original velocity field (shown here for the case of 
NGC 5457 in a plot of <v>/<x> vs. aperture number in Figure 2). While this 
may indicate a Umitation of the method, it may alternatively indicate that spiral 
patterns are in fact close to material patterns. With M. Merrifield, we are pursu- 
ing the latter possibility, while with R Rautiainen we are studying applications 
of the TW method on simulated spirals with varying pattern speeds and spiral 
strengths in order to assess the reliability of the method on spiral patterns. 
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